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Abstract 



Using a model of Frenkel excitons coupled to a bath of acoustic phonons in the 
Ch \ host medium, we study the temperature dependence of the dephasing rates and 
homogeneous line width in linear molecular aggregates. The model includes local- 
ization by disorder and predicts a power-law thermal scaling of the effective homoge- 
neous line width. The theory gives excellent agreement with temperature dependent 



t;:j- . absorption and hole-burning experiments on aggregates of the dye pseudoisocyanine. 
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1 Introduction 

Molecular J aggregates have fascinating optical properties caused by the Frenkel exciton states 
arising from the strong excitation transfer interactions in these systems. At low temperature, 
these states give rise to narrow absorption lines (exchange narrowing) [1] and ultra-fast col- 
lective spontaneous emission (exciton superradiance) [2]. The temperature dependence of the 
optical observables of molecular aggregates shows many interesting features as well, and often 
differs strongly from those of single molecules. For example, the fluorescence lifetime typically 
increases with growing temperature [2,3,4], while the fluorescence Stokes shift may exhibit an 
anomalous (non monotonous) temperature dependence by showing a growth at low tempera- 
tures [5,6]. Although the basic physics of these behaviors is often well-understood, a consistent 
quantitative explanation is mostly missing, due to the fact that the optical dynamics in these 
systems results from the complicated interplay between scattering of the excitons on disor- 
der (localization) as well as vibrations (dephasing). This also holds for the explanation of the 
experimentally observed thermal growth of the absorption line width and the homogeneous 
line width as measured in photon echo and hole burning experiments. Several authors have 
interpreted their measurements of this broadening as an activated process caused by scatter- 
ing on optical vibrations of the aggregate [3,7,8], while others concluded a power law thermal 
broadening, possibly arising from excitons scattering on vibrations in the host medium [9]. To 
add to the uncertainty, the surmised dephasing mechanisms have not been correlated with the 
temperature dependence of the fluorescence properties. 

In this paper, we model the temperature dependent dephasing rates of excitons in linear chains 
and use the results to analyze the absorption line width and hole width as a function of tem- 
perature. Our model includes scattering of the excitons on static disorder as well as on acoustic 
vibrations in the host medium. The results reveal a good agreement with experiments on ag- 
gregates of the dye pseudoisocyanine (PIC). We also comment on the successful application of 
the model to the fluorescence properties of these aggregates. 

The outline of this paper is as follows. In section 2, we describe the model. The calculation of 



the scattering rates and dephasing rates are addressed in section 3, where also the calculation 
of the absorption spectrum and the hole burning spectra are discussed. In section 4, we present 
and discuss our numerical results, while in section 5 we conclude. 

2 Model 

Our model consists of a linear Frenkel chain of A^ two- level molecules (n = 1, . . . , A^) with 
parallel transition dipoles. The Hamiltonian of this system reads 

N N 

H = J2 ^n\n){n\+ J2 ■Jnm\n){m\ , (1) 

n=l n,m=l 

where \n) denotes the state in which molecule n is excited and all others are in their ground state. 
The En are the transition energies of the individual chromophores, while the Jnm = —J/\n — m\^ 
denote the intermolecular excitation transfer interactions mediated by the transition dipoles 
(J > is the interaction between nearest neighbors). We account for disorder in the molec- 
ular energies, caused by random solvents shifts, by taking the e„ from uncorrelated Gaussian 
distributions with standard deviation a and average Eq. The exciton eigenstates (labeled u) 
follow from diagonalizing the N x N matrix with the e„ as diagonal elements and the Jnm as 
off-diagonal ones. They have energy Ei^, and site-amplitudes ipun, i-e., the eigenstates read 

N 

W) = Yl Vyn\n). (2) 

n=l 

The disorder leads to localization of the excitons state on segments of the chain; the optically 
dominant states, i.e., those with most oscillator strength, occur in the low-energy tail of the 
density of states (DOS), near the bottom of the disorder-free exciton band. 

The model includes on-site scattering of excitons on harmonic vibrations (phonons) in the host 
medium through the Hamiltonian [10] 

N 

V = YYVnq\n){n\aq: + h.c., (3) 

n=l q 

where a^ annihilates a phonon of quantum number (wave vector and branch index) g, with 
energy ujq (we set % = 1). Vnq indicates the strength of the linear exciton-phonon coupling. 



In most J aggregates, in particular for those of the prototypical PIC aggregates, the exciton- 
phonon interaction is rather weak. This is evident from the narrowness of the absorption line as 
well as the absence of a clear fluorescence Stokes shift. As a consequence, the exciton-phonon 
interaction may be accounted for in a perturbative way. This leads to transfer of population 
from exciton state v to state /i with a (scattering) rate Wy,y that may be obtained from Fermi's 
Golden Rule. The calculation of these rates is addressed in the next section. 



3 Scattering rates and spectra 

Applying Fermi's Golden Rule to the exciton-phonon scattering problem is straightforward. 
The result becomes particularly simple for a disordered host (appropriate in most experimen- 
tal situations), in which the Vnq are stochastic variables uncorrelated for the various sites n. 
Referring to Refs. [6,11] for details, we simply quote the result: 

TV 

W^. = J^{\u;,,\)Giu,.) Y: Vlndn ■ (4) 

n=l 

Here, u^^, = E^ — E^ and J^(uj) = ^nj^q |^gP^(^ ~ ^q), the one-phonon spectral density (Vg 
characterizes the coupling of this mode to the excitons). Furthermore, G{uj) = n{uj) if a; > 
and G{uj) = l + n{—uj) ifc<j < 0, with n(cj) = [exp(c(j/A;5T) — 1]^-*^, the mean thermal occupation 
number of a phonon mode of energy uj. 

The W^i, may be used to study the intra-band relaxation of excitons after photo-excitation by 
using them to set up a Pauli master equation for the occupation probabilities of the exciton 
states [6,12]. In addition, however, the W^i, govern the temperature dependent dephasing of 
the excitons, and hence their homogeneous line widths. Specifically, the thermal dephasing rate 
of the z/th exciton state is given by Fj^ = ^ Y^fi ^f^u and depends on temperature through the 
n{uj). The temperature behavior of the dephasing rates depends sensitively on the form of the 
spectral density. Here, we will assume that acoustic phonons in the host medium dominate the 
scattering process. This leads to J^{uj) = Wo{uj / J)^ , where the u^ scaling derives from the fact 
that the phonon density of states scales as uj"^ (Debye behavior) and the \Vq\'^ scales like u for 



acoustic phonons. The remaining factor Wq is a free parameter in the model, which together 
with the disorder strength a, may be used to fit experimental data. 

The above dephasing rates determine the absorption width of the individual exciton states. As 
a consequence, the total absorption spectrum reads: 

where Fi, = iJ2nVunY is the dimensionless oscillator strength of the z/th exciton state, 7,^ = 
'JqFu is its radiative rate (70 is the emission rate of a single chromophore), and the angular 
brackets denote averaging over the static disorder in the £„. In the next section we will report 
on numerical simulations of Eq. (5). In addition, the above dephasing rates may be used to 
calculate the shape and width of holes in hole burning experiments. To do this, one restricts 
the summation over z/ in eq. 5 to those states that fall within the spectral interval covered by 
the laser used to burn the hole. We will address examples of this application in section 4. 



4 Results 

It is useful to briefly consider the situation in the absence of disorder {a = 0) with only 
nearest-neighbor interactions, where the exciton states can be obtained analytically. The lowest 
{u = 1) state in the band then carries more than 80% of the oscillator strength and it can be 
shown analytically that its phonon-induced homogeneous width Fi then scales with temperature 
according to T^/^. This power law directly results from the u^ dependence of the spectral density 
and the cu^^'^ dependence of the integrated DOS close to the band bottom of one-dimensional 
systems. Inclusion of long-range dipolar interactions slightly modifles the power law to T^'®^, 
as follows from numerical calculations [13]. 

In the presence of disorder, the dephasing rates of individual exciton states vary considerably 
from one state to the other, due to stochastic fluctuations in the energy differences and wave 
function overlap summations that feature in Eq. 4 [14]. Especially at low temperature, the 
width of the distribution of these rates may be of the same order as their mean [13]. This 



makes it hard to identify a typical homogeneous width. Instead, it is more useful to directly 
simulate experimental observables. 
0.6 




Fig. 1. Calculated temperature-dependent width of the absorption spectrum (symbols) and corre- 
sponding fits to Eq. (6) (lines) for two values of the disorder strength: a = O.IJ (solid) and 0.4J 
(dashed). The fit parameters {a,p) take the values (1.2, 4.2) and (1.2, 4.3), respectively. Wq = 25 J 
and 70 = 1.5 x 10^^ J. 

For chains of 250-1000 molecules, we simulated the linear absorption spectrum according to 
Eq. 5 and extracted from it the full width at half maximum A(T) as a function of temperature. 
We have found that over an appreciable range of a and Wq values, A(T) follows an almost 
universal power law given by 



A{T) = A{0) + aWo{kBT/jy, 



(6) 



where a and p depend only weakly on the parameters a and Wq. As an example. Fig. 1 displays 
the calculated A(T) values (symbols) for two sets of parameters, together with the power-law 
fits. The (fit) parameters are given in the caption. More generally, we have found that for 
0.05J < a < 0.5J, p = 4.2 ± 0.2, while a = 1.2 ± 0.1. The scaling Eq. 6 implies that the 
total width may be interpreted as the sum of an inhomogeneous width A(0) and a dynamic 
contribution. We note that the power p lies rather close to the value of 3.85 obtained in the 
absence of disorder. The reason is that the states to which most of the population is scattered 
lie higher in the exciton band, where the DOS is not affected much by disorder. 



Interestingly, Renge and Wild [9] found that the total J-bandwidth A(T) of aggregates of 
PIC, with counter ions chloride as well as fluoride, over a wide temperature range (from 10 
K to 300 K) follows a power-law scaling as in Eq. (6). Although the power reported by these 
authors {p = 3.4) is smaller than the ones we obtained above, from direct comparison to the 
experimental data we have found that our model yields excellent fits to the experimental data 
over the entire temperature range, both for the shape and the width of the J-band. This also 
holds for the J-bandwidth of PIC-Br (bromide), measured between 1.5 K and 180 K [3]. Details 
will be published elsewhere [15]. 
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Fig. 2. Hole width F as a function of temperature measured for aggregates of PIC-I after pumping in 
the center of the J-band (triangles) [7]. The solid line is our fit; model parameters are discussed in the 
text. 



We have also calculated hole-burning spectra, according to the method described at the end of 
section 3. We present here a fit to the results reported by Hirschmann and Friedrich [7]. Using 
the hole-burning technique, they measured the homogeneous width of the exciton states in the 
center of the J-band for PIC-I (iodide) over the temperature range 350 mK to 80 K. Their data 
for the hole width F are reproduced as triangles in Fig. 2. The solid line shows our fit to these 
data, obtained by simulating disordered chains of A^ = 250 molecules. The resonant interaction 
strength and the monomer radiative rate were chosen at the accepted values of J = 600 cm^^ 
and 7o = 1.5 x lO^^J = 2.7 x 10^ s^^, respectively. Thus, the only free parameters were the 



disorder strength o and scattering strength Wq. First, we fixed the value of o by fitting the low- 
temperature (4 K) absorption spectrum, where the homogeneous broadening may be neglected. 
This yielded a = 0.21J. Next, Wq was adjusted such that the measured growth of the hole 
width was reproduced in an optimal way. Thus, we found Wq = 180J. 

We observe from Fig. 2 that our model yields a good fit to the measurements. We have also 
tried to fit the absorption and hole-burning data using a spectral density that (on average) 
scales with u according to a power different from 3 and found that is not possible to obtain a 
good fit. This suggests that scattering on acoustic phonons of the host dominates the dynamics 
of excitons in PIC aggregates. This conclusion differs markedly from previous interpretations, 
where the optical vibrations of the aggregate itself were considered as the dominant scattering 
agents, leading to sums of activated temperature dependencies [7,8]. In order to describe the 
strong increase of the line width for large temperatures, such a fit requires both a high activation 
energy (several hundreds of cm~^) and a large pre-exponential factor. The pre-factor is directly 
related to the exciton-vibration coupling constant. A simple perturbative treatment shows that 
the thus estimated coupling constant should be so large (few thousand cm^^'s) [13] that a strong 
polaron effect is to be expected. This is inconsistent with the generally accepted excitonic nature 
of the optical excitations in cyanine aggregates. 

Using the same parameters as for the for of the PIC-I hole-burning data in Fig. 2, we have 
also calculated the temperature dependence of the hole width after burning in the red and blue 
wings of the absorption line, respectively. Specifically, we used burning frequencies taken A(T) 
to the red of the peak of the absorption band and to the blue of the peak. The results are given 
in Fig. 3 together with the results when burning at the peak position. It is clearly seen that 
the hole width increases with increasing burning frequency. At low temperatures, this is to be 
expected: the dephasing of the excitons in the red wing (curve a in Fig. 3) is determined solely by 
radiative decay, whereas for the states at higher energy (curves b and c in Fig. 3) also relaxation 
to lower-energy exciton states contributes to the dephasing rate. The latter contributes grows 
with growing exciton energy. With increasing temperature, the hole width increases for all three 
burning frequencies, because downward as well as upward scattering on phonons becomes more 
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Fig. 3. Theoretical hole widths F as a function of temperature for burning frequencies in the red wing 
(dashed, a), the center (solid, b) and the blue wing (dotted, c) of the absorption spectrum. In the 
inset the arrows indicate the burning frequencies in the absorption spectrum of the aggregate at a 
temperature of 30 K. Model parameters as in Fig. 2. 

prominent. It turns out that for all cases, the hole width scales with temperature according to 
a power-like Eq. (6), where the power p depends on the burning frequency (4.1 for the red, 3.8 
for the peak, and 3.3 for the blue side). At higher temperatures, the hole widths for the three 
burning frequencies tend towards each other, because upward scattering in the exciton band 
then dominates the dephasing rates. We note that the frequency dependence of the dephasing 
rate within J-band has also been observed using photon-echo experiments, see, e.g., Ref. [17]. 



5 Concluding remarks 



In summary, we have studied the temperature dependent dephasing rate of exciton states in 
disordered Frenkel chains, caused by scattering on acoustic phonons that are characterized by a 
Debye-like (oc u^) spectral density. We have found a power-law temperature dependence for the 
absorption line width as well as the hole width obtained through hole burning. The model yields 
good quantitative explanations to absorption and hole-burning experiments on aggregates of 
the dye pseudoisocyanine with various counter-ions. 



We finally note that our model does not account for two-phonon scattering processes; in par- 
ticular we neglected pure depliasing. We have investigated such processes and found that for 
acoustic phonons their contribution to the dynamic line width of the excitons scales according 
to T'*, with s in the range 7-8.5 [13], which does not agree with the experimental findings. 
Apparently the contribution of such processes is not noticeable in aggregates of PIC. This may 
partly result from the fact that the pure-dephasing rate is suppressed by a factor that is equal 
to the delocalization size of the exciton states [13], i.e., typically by a factor of 50. 
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